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Abstract 

Semiclassical approximations to quantum dynamics are almost 
as old as quantum mechanics itself. In the approach pioneered by 
Wigner, the evolution of his quasiprobability density function on 
phase space is expressed as an asymptotic series in increasing pow- 
ers of Planck's constant, with the classical Liouvillean evolution as 
leading term. Successive semiclassical approximations to quantum 
dynamics are defined by successive terms in the series. We consider 
a complementary approach, which explores the quantum-clssical in- 
terface from the other direction. Classical dynamics is formulated 
in Hilbert space, with the Groenewold quasidensity operator as the 
image of the Liouville density on phase space. The evolution of 
the Groenewold operator is then expressed as an asymptotic series 
in increasing powers of Planck's constant. Successive semiquan- 
tum approximations to classical dynamics are defined by successive 
terms in this series, with the familiar quantum evolution as leading 
term. 
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1 INTRODUCTION 



1 Introduction 

Wigner introduced his famous quasiprobability density function on phase 
space in order to consider semiclassical approximations to the quantum 
evolution of the density matrix [1]. For a system with one degree of 
freedom and classical Hamiltonian 

H = p 2 /2m + V(q) , (1) 

he found for the evolution of the Wigner function W(q, p, t), 

W(q, p, t) t = V'(q) W(q, p, t) p - *- W(q, p, t) q 

m 



h 2 f 3v 
■— I V"'(q)W(q, p, t) ppp - JL V"(q)W(q, p, t) qpp 



+ ^-V , (q)W(q,p,t) qqp \+0(h 4 ). (2) 
m J 

Here we have introduced the notation V'(q) for dV(q)/dq and W p for 
dW/dp, etc. The first two terms on the RHS in (J2J define the classical 
Liouvillean evolution, while the terms of order h 2 define the first semiclas- 
sical approximation to the full quantum evolution of the Wigner function, 
and so on. 

Following the subsequent works of Groenewold [2] and Moyal [3] , we 
now recognize the RHS of (J2J) as the expansion in ascending powers of 
h of the star, or Moyal, bracket of the Hamiltonian H and the Wigner 
function W, 

2 3 L B R f) L 3 R 

{ H,W h = -HsMnJ/2)W, J = -----. 

h 2 

W t = {H,W} i , = HJW~—HJ 3 W + 0{h 4 ). (3) 

Here the superscripts R and L in the Janus operator J indicate the di- 
rections in which the differential operators act. The leading term H JW 
in the last equation represents the Poisson bracket of H and W, and 
corresponds to the first two terms on the RHS in Q. 

All this is very well known. It is a central ingredient of the so-called 
phase space formulation of quantum mechanics [4,5], where operators on 
Hilbert space are mapped into functions on phase space, and in particular 
the density operator is mapped into the Wigner function. 

Less well known is that, in a completely analogous way, classical me- 
chanics can be reformulated in Hilbert space [2,6,7], with the classical 
Liouville density mapped into a quasidensity operator [8-10] that we have 
called elsewhere [11] the Groenewold operator. The evolution of this op- 
erator in time is defined by what we have called [7] the odot bracket, and 



the expansion of this bracket in ascending powers of h defines a series 
of semiquantum approximations to classical dynamics, starting with the 
quantum commutator. 

In this way, we explore the classical-quantum interface in a new way, 
approaching classical mechanics from quantum mechanics, which is now 
regarded as a first approximation. So we stand on its head the traditional 
approach, which approaches quantum mechanics from classical mechanics, 
regarded as a first approximation. 



2 The Weyl-Wigner transform 

In order to see how this works, we recall firstly [4,5,12] that the phase 
space formulation of quantum mechanics is defined by the Weyl-Wigner 
transform W, which maps operators A on Hilbert space into functions A 
on phase space, 

A = W(A) , A = A(q, p) , (4) 

and in particular defines the Wigner function W — W(p/2irh) : where p is 
the density matrix defining the state of the quantum system. Then 

(A) (t) = Tv(p(t)A) = J A{q, p) W(q, p, t) dq dp = (A) (t) , 

Jw(q,p,t)dqdp=l, (5) 

but W is not in general everywhere nonnegative; it is a quasiprobability 
density function. 

In more detail, A is defined by first regarding A as an integral operator 
with kernel Ak {x, y) = {x\A\y) in the coordinate representation, and then 
setting 

A(q, p) = W(A)(q, p) = J A K (q - x/2, q + x/2) e^ x ' h dx . (6) 

The transform of the operator product on Hilbert space then defines the 
noncommutative star product on phase space, 

W(AB)=A*B, (7) 

leading to the star bracket as the image of (1/ihx) the commutator, 

W([i, B]/ih) = (A*B-B* A)/ih = {A, B}± . (8) 



It can now be seen that the quantum evolution of the density matrix 
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3 THE ODOT PRODUCT AND ODOT BRACKET 



is mapped by the Weyl-Wigner transform W into the evolution equation 
for the Wigner function 

W t = {H, W}± (10) 

as in 10, so leading to the sequence of semiclassical approximations as 
described in the Introduction. In if is the quantum Hamiltonian 
operator, so that H = W(H). 

In order to define semiquantum approximations to classical dynamics, 
we begin by considering the inverse Weyl-Wigner transform W _1 , which 
maps functions A on phase space into operators A on Hilbert space, so 
enabling a Hilbert space formulation of classical mechanics [6,7]. We have 

A K ( X ,y)=W- 1 (A) K (x,y) = 7 ±- f A([x + y}/2, p) e *(*-»>/» dp , (11) 

which defines the kernel of A, and hence A itself, in terms of A. In 
particular the Groenewold operator G(t) is defined as G — W~ 1 (2irhp). 
Then 

(A)(t) = J A{q, p)dqdp = Tr(AG(t)) = (A) , 

Tr(G(t))=l, (12) 

but G is not nonnegative definite in general; it is a quasidensity operator. 

It can be seen that the development so far is completely analogous to 
the development of the phase space formulation of quantum mechanics. 
We can use W _1 to map all of classical mechanics into a Hilbert space 
formulation [6,7]. To complete the story, we need to say what happens to 
the classical evolution of the Liouville density 

Pt = HJp (13) 

under the action of W _1 . Obviously the LHS maps into Gt\ the question 
is, what happens to the Poisson bracket on the RHS. 

3 The odot product and odot bracket 

To proceed, we note firstly by analogy with the definition of the star 
product that we can define a commutative odot product of operators on 
Hilbert space 

A(DB =W- 1 (AB). (14) 

This is an interesting product, quite distinct from the well known Jor- 
dan product of operators, which is also commutative. Unlike the Jordan 
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product, however, this odot product is associative. Some of its other 
characteristic properties have been described elsewhere [6,7,13]. 
Next we note that A q — {A, p}+ so that 

W-\A q ) = ^[p 1 A]=A qi say, (15) 

where p = W~ 1 (p). Similarly, we define 

W-\A P ) = q]=A p , 

W-\A qp ) = -^[[p, A],q]=A qp , etc. (16) 

Then 

W- 1 (AJB)=A q QB p -A p QB q = ^-[A,B] e , say, (17) 

which defines the odot bracket as the inverse image of ihx the Poisson 
bracket. From the Poisson bracket it inherits antisymmetry and a Jacobi 
identity. 

Our next task is to find an expansion of the odot bracket in ascending 
powers of K, analogous to the expansion J2J of the star bracket, in order to 
define a sequence of semiquantum approximations to classical dynamics. 



4 Semiquantum mechanics 

We set 

M = fsin(f), (18) 

so we can write for any two functions A and B, 

{A, B}+ = AM B , (19) 

and note therefore that 

W'UAM B) = ^-[A, B}. (20) 
in 

Next we write the Poisson bracket as 

A3 B = A hJ/2 , - MB (21) 

sin(ftJ/2) V ; 

and, noting that 

= 1 + 6 2 /6 + 76» 4 /360 - . . . (22) 



sin(0) 



G 



5 EXAMPLES 



we obtain from i|2U|) and lltifl that 
W^iAJB) = 

±[A, B}- 1 ^ ([A qq , B pp ] - 2[A qp , B qp ] + [A pp , B qq }) + 0(ft 3 ) . (23) 

Now we can answer the question raised at the end of Section 2. Applying 
to both sides of Ijl3|l . we obtain 

Gt = 

G] - ^ ([&„, G pp ] ~ 2[H qp , G qp ] + [H pp , G qq \) + 0(ft 3 ) . (24) 

Thus we see the evolution of G in Hilbert space, which is equivalent to 
the classical evolution of the Liouville density p in phase space, is given as 
a series in ascending powers of ft. Keeping successively more terms in the 
series, we define a sequence of approximations to the classical evolution. 
Note that the lowest order term is just the quantum evolution, which now 
appears as the lowest order approximation to classical dynamics. As we 
add more and more terms, we have the possibility to explore the classical- 
quantum interface starting from the quantum side. This is completely 
complementary to what we normally do with semiclassical approximations 
to quantum mechanics. 

5 Examples 

J.G. Wood and I have explored semiquantum and semiclassical approx- 
imations for simple nonlinear systems with one degree of freedom [14]. 
Rather than Hamiltonians of the form Q , we considered 

K 

H = Ej2h(H /E) k , (25) 

fc=0 

where Hq is the simple harmonic oscillator Hamiltonian 

H a = p 2 /2m + muj 2 q 2 /2 , (26) 

and E, bk are constants. These have the advantage that they are an- 
alytically tractable, but still show characteristic differences between the 
classical and quantum evolutions [15]. In particular we considered 

H 2 = Hl/E , H 2 = [ifuj {N 2 + N + 1/4) , 

H 3 = Hl/E 2 , H 3 = p 2 hw {N 3 + 3N 2 /2 + 2N + 3/4) , 



where H = hw (N + 1 /2) , /j,= tw/E, (27) 
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with H = W _1 (ii") in each case, and N the usual oscillator number 
operator. 

As an initial Liouville density on phase space, we took a Gaussian 
p for which the initial Groenewold operator G = W~ 1 (2Trhp) equals a 
true density operator, namely the density operator for a pure coherent 
state. Differences between the classical and quantum evolutions of such 
an initial state, with the Hamiltonians if 2 and H2, respectively, are im- 
mediately apparent in the phase space plots Fig. ^ an d Fig. Under 
the classical evolution, the density stays positive everywhere, but devel- 
ops "whorls," whereas under the quantum evolution, the density (Wigner 
function) becomes negative on some regions (shown in white) and is pe- 
riodic [15]. Conversely, under the classical dynamics, the Groenewold 
operator G develops negative eigenvalues [8-11,16], whereas under the 
quantum evolution, such an initial pure-state density operator stays pos- 
itive definite, with eigenvalues and 1. Simliar remarks apply with the 
Hamiltonians H3 and H3. 

There is no nontrivial semiquantum or semiclassical approximation 
"in between" the classical and quantum dynamics for the Hamiltonians 
H2 and H2. Each of the series J3J) and (|24J) has just two terms in this 
case. For the series (0), the leading term defines the classical evolution, 
and adding the next term produces the full quantum dynamics. Similarly, 
for the series Il24fl . the leading term defines the quantum evolution, and 
adding the next term produces the full classical dynamics. 

To see interesting differences betweem semiclassical and semiquantum 
approximations, we considered H% and H3, for which there are three terms 
in each of the series (01 and (|24|l . Thus we can compare the classical 
evolution and the first semiclassical approximation to quantum dynamics, 
and also the quantum evolution and the first semiquantum approximation 
to classical dynamics. In Fig. |3|we plot the expectation values of q, p in the 
classical and semiclassical cases, calculated at each time using the Liouville 
density p or the Wigner function W as in ©, and also the expectation 
values of q and p in the quantum and semiquantum cases, calculated 
at each time using the density operator p or the Groenewold operator 
G as in 112fl . From our results it is already clear that semiquantum 
and semiclassical approximations provide different information about the 
interface between classical and quantum behaviours. 

These expectation values compare "classical-like" properties of the dif- 
ferent evolutions. We also considered "quantum- like" properties of the 
diferent evolutions, in particular the largest and smallest eigenvalues of 
W^ 1 (p) and W~ 1 (W) in the classical and semiclassical cases, and the 
largest and smallest eigenvalues of p and G in the quantum and semi- 
quantum cases. The results are shown in Fig. 
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5 EXAMPLES 




Figure 1: Density plots showing the classical evolution of an initial Gaus- 
sian density centered at qo — 0.5, po ~ as generated by the Hamiltonian 
H = Hq/E. The parameters m,uj,E have been set equal to 1, and the 
times of the plots are, from left to right and top to bottom, t — ir/4, 
t = 7r/2, t = 3tt/4 and t = it. 
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Figure 2: Quantum evolution of an initial Gaussian Wigner function, with 
the same parameter values used in Figure^ an d shown at the same times. 
Regions on which the pseudo-density becomes negative are shown in white 
in the first three plots. 
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Figure 3: Comparison of first moments of a = (q + ip)/y2 for classi- 
cal, semiquantum, quantum and semiclassical evolutions generated by the 
Hamiltonian H = Hq/E 2 . Points on the classical, quantum, semiclassical 
and semiquantum curves are labelled by +, x, o and □, respectively. The 
evolution is over the time- interval [0, 7r] and again m = u> = E = 1, with 
H = 1/2, a = 0.5. 
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Figure 4: Comparison of largest two and least two eigenvalues for, from 
left to right and top to bottom, classical, semiquantum and semiclassical 
evolutions generated by H = Hq/E 2 , for the time interval [0, it], and with 
m = lu = E = 1 and /i = 1/2, ao = 0.5. The evolution of the first moment 
is reproduced from Fig. [3] in the graph at bottom right for comparison. 
Each of the other graphs also features the quantum spectrum {0, 1} and in 
all graphs the values at the time-points t = 1, 2, 3 are marked (classical), 
o (semiclassical), A (quantum) and □ (semiquantum). 
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6 Conclusions 

Semiquantum mechanics opens a new window on the interface between 
classical and quantum mechanics. Our investigations of examples as de- 
scribed above have not yet gone very far, but already we can say that 
semiquantum approximations show characteristic differences from semi- 
classical approximations for a given nonlinear system. More details are 
given in Ref. [14]. 

It is important to explore the nature of semiquantum approximations 
for other systems. For example, from the classical Hamiltonian H as in 
JU, wc obtain 

H = f/2m + V(q) , (28) 

and substituting in l|24|) . we get for the evolution of the Groenewold op- 
erator in such cases 

G t = ^[f/2m + V(q), G] - ^[V"(q), G pp ] + 0(h 3 ) , (29) 

which is to be compared with Wigner's famous formula JSJ for the evolu- 
tion of the Wigner function. The implications of l|29|l for various impor- 
tant choices of V, in particular exactly solvable cases, should be examined. 

Of even more interest of course are systems with more degrees of free- 
dom that show chaos at the classical level. For such systems it should be 
particularly interesting to consider differences in the spectral properties 
of G at different times in different approximations, in classically chaotic 
or integrable regimes, as well as the different behaviours of expectation 
values of important phase space variables. 

We hope to investigate some of these problems. 

Acknowledgments 

Thanks to Dr. James Wood for much advice and assistance, and to a 
referee for bringing the paper by Vercin [6] to our attention. This work 
was supported by Australian Research Council Grant DP0450778. 

References 

[1] E.P. Wigner, Phys. Rev. 40, 749 (1932). 

[2] H. Groenewold, Physica 12, 406 (1946). 

[3] J. Moyal, Proc. Camb. Phil. Soc. 45, 99 (1949). 

[4] D.A. Dubin, M.A. Hennings and T.B. Smith, Mathematical Aspects 
of Weyl Quantization and Phase (Singapore, World Scientific, 2000). 



REFERENCES 



13 



[5] C. Zachos, Int. J. Mod. Phys. 17, 297 (2002). 

[6] A. Vercin, Int. J. Theoret. Phys. 39, 2063 (2000). 

[7] A.J. Bracken, J. Phys. A: Math. Gen. 36, L329 (2003). 

[8] J.G. Muga and R.F. Snider, Europhys. Lett. 19, 569 (1992). 

[9] J.G. Muga, R. Sala, and R.F. Snider, Physica Scripta 47, 732 (1993). 

[10] R. Sala and J.G. Muga, Phys. Lett. A 192, 180 (1994). 

[11] A.J. Bracken and J.G. Wood, Europhys. Lett. 68, 1 (2004). 

[12] A.J. Bracken, G. Cassinelli and J.G. Wood, J. Phys. A: Math. Gen. 
36, 1033 (2003). 

[13] D.A. Dubin and M.A. Hennings, J. Phys. A: Math. Gen. 37, 6693 
(2004). 

[14] A.J. Bracken and J.G. Wood, Semiquantum versus semiclassical 
mechanics for simple nonlinear systems, University of Queensland 
preprint, quant-ph/0511227 (to appear in Phys. Rev. A.) 

[15] G.J. Milburn, Phys. Rev. A 33, 674 (1986). 

[16] S. Habib, K. Jacobs, H. Mabuchi, R. Ryne, K. Shizume and B. Sun- 
daram, Phys. Rev. Lett. 88, 040402 (2002). 



